3.4.82 \(\int \frac {\sqrt {d+e x}}{(b x+c x^2)^3} \, dx\) [382]

Optimal. Leaf size=245 \[ -\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (b (c d-b e) (12 c d-b e)+c \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x\right )}{4 b^4 d (c d-b e) \left (b x+c x^2\right )}-\frac {\left (48 c^2 d^2-12 b c d e-b^2 e^2\right ) \tanh ^{-1}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{4 b^5 d^{3/2}}+\frac {c^{3/2} \left (48 c^2 d^2-84 b c d e+35 b^2 e^2\right ) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{4 b^5 (c d-b e)^{3/2}} \]

[Out]

-1/4*(-b^2*e^2-12*b*c*d*e+48*c^2*d^2)*arctanh((e*x+d)^(1/2)/d^(1/2))/b^5/d^(3/2)+1/4*c^(3/2)*(35*b^2*e^2-84*b*
c*d*e+48*c^2*d^2)*arctanh(c^(1/2)*(e*x+d)^(1/2)/(-b*e+c*d)^(1/2))/b^5/(-b*e+c*d)^(3/2)-1/2*(2*c*x+b)*(e*x+d)^(
1/2)/b^2/(c*x^2+b*x)^2+1/4*(b*(-b*e+c*d)*(-b*e+12*c*d)+c*(b^2*e^2-24*b*c*d*e+24*c^2*d^2)*x)*(e*x+d)^(1/2)/b^4/
d/(-b*e+c*d)/(c*x^2+b*x)

________________________________________________________________________________________

Rubi [A]
time = 0.27, antiderivative size = 245, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, integrand size = 21, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.238, Rules used = {750, 836, 840, 1180, 214} \begin {gather*} -\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}-\frac {\left (-b^2 e^2-12 b c d e+48 c^2 d^2\right ) \tanh ^{-1}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{4 b^5 d^{3/2}}+\frac {c^{3/2} \left (35 b^2 e^2-84 b c d e+48 c^2 d^2\right ) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{4 b^5 (c d-b e)^{3/2}}+\frac {\sqrt {d+e x} \left (c x \left (b^2 e^2-24 b c d e+24 c^2 d^2\right )+b (c d-b e) (12 c d-b e)\right )}{4 b^4 d \left (b x+c x^2\right ) (c d-b e)} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[Sqrt[d + e*x]/(b*x + c*x^2)^3,x]

[Out]

-1/2*((b + 2*c*x)*Sqrt[d + e*x])/(b^2*(b*x + c*x^2)^2) + (Sqrt[d + e*x]*(b*(c*d - b*e)*(12*c*d - b*e) + c*(24*
c^2*d^2 - 24*b*c*d*e + b^2*e^2)*x))/(4*b^4*d*(c*d - b*e)*(b*x + c*x^2)) - ((48*c^2*d^2 - 12*b*c*d*e - b^2*e^2)
*ArcTanh[Sqrt[d + e*x]/Sqrt[d]])/(4*b^5*d^(3/2)) + (c^(3/2)*(48*c^2*d^2 - 84*b*c*d*e + 35*b^2*e^2)*ArcTanh[(Sq
rt[c]*Sqrt[d + e*x])/Sqrt[c*d - b*e]])/(4*b^5*(c*d - b*e)^(3/2))

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 750

Int[((d_.) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(d + e*x)^m*(b + 2*c
*x)*((a + b*x + c*x^2)^(p + 1)/((p + 1)*(b^2 - 4*a*c))), x] - Dist[1/((p + 1)*(b^2 - 4*a*c)), Int[(d + e*x)^(m
- 1)*(b*e*m + 2*c*d*(2*p + 3) + 2*c*e*(m + 2*p + 3)*x)*(a + b*x + c*x^2)^(p + 1), x], x] /; FreeQ[{a, b, c, d
, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && NeQ[2*c*d - b*e, 0] && LtQ[p, -1] && GtQ[m
, 0] && (LtQ[m, 1] || (ILtQ[m + 2*p + 3, 0] && NeQ[m, 2])) && IntQuadraticQ[a, b, c, d, e, m, p, x]

Rule 836

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp
[(d + e*x)^(m + 1)*(f*(b*c*d - b^2*e + 2*a*c*e) - a*g*(2*c*d - b*e) + c*(f*(2*c*d - b*e) - g*(b*d - 2*a*e))*x)
*((a + b*x + c*x^2)^(p + 1)/((p + 1)*(b^2 - 4*a*c)*(c*d^2 - b*d*e + a*e^2))), x] + Dist[1/((p + 1)*(b^2 - 4*a*
c)*(c*d^2 - b*d*e + a*e^2)), Int[(d + e*x)^m*(a + b*x + c*x^2)^(p + 1)*Simp[f*(b*c*d*e*(2*p - m + 2) + b^2*e^2
*(p + m + 2) - 2*c^2*d^2*(2*p + 3) - 2*a*c*e^2*(m + 2*p + 3)) - g*(a*e*(b*e - 2*c*d*m + b*e*m) - b*d*(3*c*d -
b*e + 2*c*d*p - b*e*p)) + c*e*(g*(b*d - 2*a*e) - f*(2*c*d - b*e))*(m + 2*p + 4)*x, x], x], x] /; FreeQ[{a, b,
c, d, e, f, g, m}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && LtQ[p, -1] && (IntegerQ[m] ||
IntegerQ[p] || IntegersQ[2*m, 2*p])

Rule 840

Int[((f_.) + (g_.)*(x_))/(Sqrt[(d_.) + (e_.)*(x_)]*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)), x_Symbol] :> Dist[2,
Subst[Int[(e*f - d*g + g*x^2)/(c*d^2 - b*d*e + a*e^2 - (2*c*d - b*e)*x^2 + c*x^4), x], x, Sqrt[d + e*x]], x] /
; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0]

Rule 1180

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rubi steps

\begin {align*} \int \frac {\sqrt {d+e x}}{\left (b x+c x^2\right )^3} \, dx &=-\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}+\frac {\int \frac {-6 c d+\frac {b e}{2}-5 c e x}{\sqrt {d+e x} \left (b x+c x^2\right )^2} \, dx}{2 b^2}\\ &=-\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (b (c d-b e) (12 c d-b e)+c \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x\right )}{4 b^4 d (c d-b e) \left (b x+c x^2\right )}-\frac {\int \frac {-\frac {1}{4} (c d-b e) \left (48 c^2 d^2-12 b c d e-b^2 e^2\right )-\frac {1}{4} c e \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x}{\sqrt {d+e x} \left (b x+c x^2\right )} \, dx}{2 b^4 d (c d-b e)}\\ &=-\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (b (c d-b e) (12 c d-b e)+c \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x\right )}{4 b^4 d (c d-b e) \left (b x+c x^2\right )}-\frac {\text {Subst}\left (\int \frac {-\frac {1}{4} e (c d-b e) \left (48 c^2 d^2-12 b c d e-b^2 e^2\right )+\frac {1}{4} c d e \left (24 c^2 d^2-24 b c d e+b^2 e^2\right )-\frac {1}{4} c e \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x^2}{c d^2-b d e+(-2 c d+b e) x^2+c x^4} \, dx,x,\sqrt {d+e x}\right )}{b^4 d (c d-b e)}\\ &=-\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (b (c d-b e) (12 c d-b e)+c \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x\right )}{4 b^4 d (c d-b e) \left (b x+c x^2\right )}+\frac {\left (c \left (48 c^2 d^2-12 b c d e-b^2 e^2\right )\right ) \text {Subst}\left (\int \frac {1}{-\frac {b e}{2}+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{4 b^5 d}-\frac {\left (c^2 \left (48 c^2 d^2-84 b c d e+35 b^2 e^2\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b e}{2}+\frac {1}{2} (-2 c d+b e)+c x^2} \, dx,x,\sqrt {d+e x}\right )}{4 b^5 (c d-b e)}\\ &=-\frac {(b+2 c x) \sqrt {d+e x}}{2 b^2 \left (b x+c x^2\right )^2}+\frac {\sqrt {d+e x} \left (b (c d-b e) (12 c d-b e)+c \left (24 c^2 d^2-24 b c d e+b^2 e^2\right ) x\right )}{4 b^4 d (c d-b e) \left (b x+c x^2\right )}-\frac {\left (48 c^2 d^2-12 b c d e-b^2 e^2\right ) \tanh ^{-1}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{4 b^5 d^{3/2}}+\frac {c^{3/2} \left (48 c^2 d^2-84 b c d e+35 b^2 e^2\right ) \tanh ^{-1}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {c d-b e}}\right )}{4 b^5 (c d-b e)^{3/2}}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]
time = 1.36, size = 247, normalized size = 1.01 \begin {gather*} \frac {\frac {b \sqrt {d+e x} \left (24 c^4 d^2 x^3+12 b c^3 d x^2 (3 d-2 e x)+b^4 e (2 d+e x)+b^2 c^2 x \left (8 d^2-37 d e x+e^2 x^2\right )+b^3 c \left (-2 d^2-9 d e x+2 e^2 x^2\right )\right )}{d (c d-b e) x^2 (b+c x)^2}+\frac {c^{3/2} \left (48 c^2 d^2-84 b c d e+35 b^2 e^2\right ) \tan ^{-1}\left (\frac {\sqrt {c} \sqrt {d+e x}}{\sqrt {-c d+b e}}\right )}{(-c d+b e)^{3/2}}+\frac {\left (-48 c^2 d^2+12 b c d e+b^2 e^2\right ) \tanh ^{-1}\left (\frac {\sqrt {d+e x}}{\sqrt {d}}\right )}{d^{3/2}}}{4 b^5} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[Sqrt[d + e*x]/(b*x + c*x^2)^3,x]

[Out]

((b*Sqrt[d + e*x]*(24*c^4*d^2*x^3 + 12*b*c^3*d*x^2*(3*d - 2*e*x) + b^4*e*(2*d + e*x) + b^2*c^2*x*(8*d^2 - 37*d
*e*x + e^2*x^2) + b^3*c*(-2*d^2 - 9*d*e*x + 2*e^2*x^2)))/(d*(c*d - b*e)*x^2*(b + c*x)^2) + (c^(3/2)*(48*c^2*d^
2 - 84*b*c*d*e + 35*b^2*e^2)*ArcTan[(Sqrt[c]*Sqrt[d + e*x])/Sqrt[-(c*d) + b*e]])/(-(c*d) + b*e)^(3/2) + ((-48*
c^2*d^2 + 12*b*c*d*e + b^2*e^2)*ArcTanh[Sqrt[d + e*x]/Sqrt[d]])/d^(3/2))/(4*b^5)

________________________________________________________________________________________

Maple [A]
time = 0.50, size = 258, normalized size = 1.05

method	result	size
derivativedivides	\(2 e^{5} \left (-\frac {\frac {\frac {b e \left (b e -12 c d \right ) \left (e x +d \right )^{\frac {3}{2}}}{8 d}+\left (\frac {1}{8} b^{2} e^{2}+\frac {3}{2} b c d e \right ) \sqrt {e x +d}}{e^{2} x^{2}}-\frac {\left (b^{2} e^{2}+12 b c d e -48 d^{2} c^{2}\right ) \arctanh \left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{8 d^{\frac {3}{2}}}}{b^{5} e^{5}}+\frac {c^{2} \left (\frac {\frac {c b e \left (11 b e -12 c d \right ) \left (e x +d \right )^{\frac {3}{2}}}{8 b e -8 c d}+\frac {b e \left (13 b e -12 c d \right ) \sqrt {e x +d}}{8}}{\left (c \left (e x +d \right )+b e -c d \right )^{2}}+\frac {\left (35 b^{2} e^{2}-84 b c d e +48 d^{2} c^{2}\right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{8 \left (b e -c d \right ) \sqrt {\left (b e -c d \right ) c}}\right )}{b^{5} e^{5}}\right )\)	\(258\)
default	\(2 e^{5} \left (-\frac {\frac {\frac {b e \left (b e -12 c d \right ) \left (e x +d \right )^{\frac {3}{2}}}{8 d}+\left (\frac {1}{8} b^{2} e^{2}+\frac {3}{2} b c d e \right ) \sqrt {e x +d}}{e^{2} x^{2}}-\frac {\left (b^{2} e^{2}+12 b c d e -48 d^{2} c^{2}\right ) \arctanh \left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{8 d^{\frac {3}{2}}}}{b^{5} e^{5}}+\frac {c^{2} \left (\frac {\frac {c b e \left (11 b e -12 c d \right ) \left (e x +d \right )^{\frac {3}{2}}}{8 b e -8 c d}+\frac {b e \left (13 b e -12 c d \right ) \sqrt {e x +d}}{8}}{\left (c \left (e x +d \right )+b e -c d \right )^{2}}+\frac {\left (35 b^{2} e^{2}-84 b c d e +48 d^{2} c^{2}\right ) \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{8 \left (b e -c d \right ) \sqrt {\left (b e -c d \right ) c}}\right )}{b^{5} e^{5}}\right )\)	\(258\)
risch	\(-\frac {\sqrt {e x +d}\, \left (b e x -12 c d x +2 b d \right )}{4 d \,b^{4} x^{2}}+\frac {11 e^{2} c^{3} \left (e x +d \right )^{\frac {3}{2}}}{4 b^{3} \left (c e x +b e \right )^{2} \left (b e -c d \right )}-\frac {3 d e \,c^{4} \left (e x +d \right )^{\frac {3}{2}}}{b^{4} \left (c e x +b e \right )^{2} \left (b e -c d \right )}+\frac {13 e^{2} c^{2} \sqrt {e x +d}}{4 b^{3} \left (c e x +b e \right )^{2}}-\frac {3 d e \,c^{3} \sqrt {e x +d}}{b^{4} \left (c e x +b e \right )^{2}}+\frac {35 e^{2} c^{2} \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{4 b^{3} \left (b e -c d \right ) \sqrt {\left (b e -c d \right ) c}}-\frac {21 d e \,c^{3} \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{b^{4} \left (b e -c d \right ) \sqrt {\left (b e -c d \right ) c}}+\frac {12 d^{2} c^{4} \arctan \left (\frac {c \sqrt {e x +d}}{\sqrt {\left (b e -c d \right ) c}}\right )}{b^{5} \left (b e -c d \right ) \sqrt {\left (b e -c d \right ) c}}+\frac {e^{2} \arctanh \left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right )}{4 b^{3} d^{\frac {3}{2}}}+\frac {3 e \arctanh \left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right ) c}{b^{4} \sqrt {d}}-\frac {12 \sqrt {d}\, \arctanh \left (\frac {\sqrt {e x +d}}{\sqrt {d}}\right ) c^{2}}{b^{5}}\)	\(396\)

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((e*x+d)^(1/2)/(c*x^2+b*x)^3,x,method=_RETURNVERBOSE)

[Out]

2*e^5*(-1/b^5/e^5*((1/8*b*e*(b*e-12*c*d)/d*(e*x+d)^(3/2)+(1/8*b^2*e^2+3/2*b*c*d*e)*(e*x+d)^(1/2))/e^2/x^2-1/8*
(b^2*e^2+12*b*c*d*e-48*c^2*d^2)/d^(3/2)*arctanh((e*x+d)^(1/2)/d^(1/2)))+1/b^5/e^5*c^2*((1/8*c*b*e*(11*b*e-12*c
*d)/(b*e-c*d)*(e*x+d)^(3/2)+1/8*b*e*(13*b*e-12*c*d)*(e*x+d)^(1/2))/(c*(e*x+d)+b*e-c*d)^2+1/8*(35*b^2*e^2-84*b*
c*d*e+48*c^2*d^2)/(b*e-c*d)/((b*e-c*d)*c)^(1/2)*arctan(c*(e*x+d)^(1/2)/((b*e-c*d)*c)^(1/2))))

________________________________________________________________________________________

Maxima [F(-2)]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Exception raised: ValueError} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)^(1/2)/(c*x^2+b*x)^3,x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(c*d-%e*b>0)', see `assume?` fo
r more detai

________________________________________________________________________________________

Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 583 vs. \(2 (226) = 452\).
time = 4.10, size = 2372, normalized size = 9.68 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)^(1/2)/(c*x^2+b*x)^3,x, algorithm="fricas")

[Out]

[-1/8*((48*c^5*d^4*x^4 + 96*b*c^4*d^4*x^3 + 48*b^2*c^3*d^4*x^2 + 35*(b^2*c^3*d^2*x^4 + 2*b^3*c^2*d^2*x^3 + b^4
*c*d^2*x^2)*e^2 - 84*(b*c^4*d^3*x^4 + 2*b^2*c^3*d^3*x^3 + b^3*c^2*d^3*x^2)*e)*sqrt(c/(c*d - b*e))*log((2*c*d -
2*(c*d - b*e)*sqrt(x*e + d)*sqrt(c/(c*d - b*e)) + (c*x - b)*e)/(c*x + b)) + (48*c^5*d^3*x^4 + 96*b*c^4*d^3*x^
3 + 48*b^2*c^3*d^3*x^2 + (b^3*c^2*x^4 + 2*b^4*c*x^3 + b^5*x^2)*e^3 + 11*(b^2*c^3*d*x^4 + 2*b^3*c^2*d*x^3 + b^4
*c*d*x^2)*e^2 - 60*(b*c^4*d^2*x^4 + 2*b^2*c^3*d^2*x^3 + b^3*c^2*d^2*x^2)*e)*sqrt(d)*log((x*e + 2*sqrt(x*e + d)
*sqrt(d) + 2*d)/x) - 2*(24*b*c^4*d^3*x^3 + 36*b^2*c^3*d^3*x^2 + 8*b^3*c^2*d^3*x - 2*b^4*c*d^3 + (b^3*c^2*d*x^3
+ 2*b^4*c*d*x^2 + b^5*d*x)*e^2 - (24*b^2*c^3*d^2*x^3 + 37*b^3*c^2*d^2*x^2 + 9*b^4*c*d^2*x - 2*b^5*d^2)*e)*sqr
t(x*e + d))/(b^5*c^3*d^3*x^4 + 2*b^6*c^2*d^3*x^3 + b^7*c*d^3*x^2 - (b^6*c^2*d^2*x^4 + 2*b^7*c*d^2*x^3 + b^8*d^
2*x^2)*e), 1/8*(2*(48*c^5*d^4*x^4 + 96*b*c^4*d^4*x^3 + 48*b^2*c^3*d^4*x^2 + 35*(b^2*c^3*d^2*x^4 + 2*b^3*c^2*d^
2*x^3 + b^4*c*d^2*x^2)*e^2 - 84*(b*c^4*d^3*x^4 + 2*b^2*c^3*d^3*x^3 + b^3*c^2*d^3*x^2)*e)*sqrt(-c/(c*d - b*e))*
arctan(-(c*d - b*e)*sqrt(x*e + d)*sqrt(-c/(c*d - b*e))/(c*x*e + c*d)) - (48*c^5*d^3*x^4 + 96*b*c^4*d^3*x^3 + 4
8*b^2*c^3*d^3*x^2 + (b^3*c^2*x^4 + 2*b^4*c*x^3 + b^5*x^2)*e^3 + 11*(b^2*c^3*d*x^4 + 2*b^3*c^2*d*x^3 + b^4*c*d*
x^2)*e^2 - 60*(b*c^4*d^2*x^4 + 2*b^2*c^3*d^2*x^3 + b^3*c^2*d^2*x^2)*e)*sqrt(d)*log((x*e + 2*sqrt(x*e + d)*sqrt
(d) + 2*d)/x) + 2*(24*b*c^4*d^3*x^3 + 36*b^2*c^3*d^3*x^2 + 8*b^3*c^2*d^3*x - 2*b^4*c*d^3 + (b^3*c^2*d*x^3 + 2*
b^4*c*d*x^2 + b^5*d*x)*e^2 - (24*b^2*c^3*d^2*x^3 + 37*b^3*c^2*d^2*x^2 + 9*b^4*c*d^2*x - 2*b^5*d^2)*e)*sqrt(x*e
+ d))/(b^5*c^3*d^3*x^4 + 2*b^6*c^2*d^3*x^3 + b^7*c*d^3*x^2 - (b^6*c^2*d^2*x^4 + 2*b^7*c*d^2*x^3 + b^8*d^2*x^2
)*e), 1/8*(2*(48*c^5*d^3*x^4 + 96*b*c^4*d^3*x^3 + 48*b^2*c^3*d^3*x^2 + (b^3*c^2*x^4 + 2*b^4*c*x^3 + b^5*x^2)*e
^3 + 11*(b^2*c^3*d*x^4 + 2*b^3*c^2*d*x^3 + b^4*c*d*x^2)*e^2 - 60*(b*c^4*d^2*x^4 + 2*b^2*c^3*d^2*x^3 + b^3*c^2*
d^2*x^2)*e)*sqrt(-d)*arctan(sqrt(x*e + d)*sqrt(-d)/d) - (48*c^5*d^4*x^4 + 96*b*c^4*d^4*x^3 + 48*b^2*c^3*d^4*x^
2 + 35*(b^2*c^3*d^2*x^4 + 2*b^3*c^2*d^2*x^3 + b^4*c*d^2*x^2)*e^2 - 84*(b*c^4*d^3*x^4 + 2*b^2*c^3*d^3*x^3 + b^3
*c^2*d^3*x^2)*e)*sqrt(c/(c*d - b*e))*log((2*c*d - 2*(c*d - b*e)*sqrt(x*e + d)*sqrt(c/(c*d - b*e)) + (c*x - b)*
e)/(c*x + b)) + 2*(24*b*c^4*d^3*x^3 + 36*b^2*c^3*d^3*x^2 + 8*b^3*c^2*d^3*x - 2*b^4*c*d^3 + (b^3*c^2*d*x^3 + 2*
b^4*c*d*x^2 + b^5*d*x)*e^2 - (24*b^2*c^3*d^2*x^3 + 37*b^3*c^2*d^2*x^2 + 9*b^4*c*d^2*x - 2*b^5*d^2)*e)*sqrt(x*e
+ d))/(b^5*c^3*d^3*x^4 + 2*b^6*c^2*d^3*x^3 + b^7*c*d^3*x^2 - (b^6*c^2*d^2*x^4 + 2*b^7*c*d^2*x^3 + b^8*d^2*x^2
)*e), 1/4*((48*c^5*d^4*x^4 + 96*b*c^4*d^4*x^3 + 48*b^2*c^3*d^4*x^2 + 35*(b^2*c^3*d^2*x^4 + 2*b^3*c^2*d^2*x^3 +
b^4*c*d^2*x^2)*e^2 - 84*(b*c^4*d^3*x^4 + 2*b^2*c^3*d^3*x^3 + b^3*c^2*d^3*x^2)*e)*sqrt(-c/(c*d - b*e))*arctan(
-(c*d - b*e)*sqrt(x*e + d)*sqrt(-c/(c*d - b*e))/(c*x*e + c*d)) + (48*c^5*d^3*x^4 + 96*b*c^4*d^3*x^3 + 48*b^2*c
^3*d^3*x^2 + (b^3*c^2*x^4 + 2*b^4*c*x^3 + b^5*x^2)*e^3 + 11*(b^2*c^3*d*x^4 + 2*b^3*c^2*d*x^3 + b^4*c*d*x^2)*e^
2 - 60*(b*c^4*d^2*x^4 + 2*b^2*c^3*d^2*x^3 + b^3*c^2*d^2*x^2)*e)*sqrt(-d)*arctan(sqrt(x*e + d)*sqrt(-d)/d) + (2
4*b*c^4*d^3*x^3 + 36*b^2*c^3*d^3*x^2 + 8*b^3*c^2*d^3*x - 2*b^4*c*d^3 + (b^3*c^2*d*x^3 + 2*b^4*c*d*x^2 + b^5*d*
x)*e^2 - (24*b^2*c^3*d^2*x^3 + 37*b^3*c^2*d^2*x^2 + 9*b^4*c*d^2*x - 2*b^5*d^2)*e)*sqrt(x*e + d))/(b^5*c^3*d^3*
x^4 + 2*b^6*c^2*d^3*x^3 + b^7*c*d^3*x^2 - (b^6*c^2*d^2*x^4 + 2*b^7*c*d^2*x^3 + b^8*d^2*x^2)*e)]

________________________________________________________________________________________

Sympy [B] Leaf count of result is larger than twice the leaf count of optimal. 2584 vs. \(2 (224) = 448\).
time = 41.12, size = 2584, normalized size = 10.55 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)**(1/2)/(c*x**2+b*x)**3,x)

[Out]

10*c**4*d**2*e**2*sqrt(d + e*x)/(8*b**7*e**4 - 16*b**6*c*d*e**3 + 16*b**6*c*e**4*x - 48*b**5*c**2*d*e**3*x + 8
*b**5*c**2*e**2*(d + e*x)**2 + 16*b**4*c**3*d**3*e + 48*b**4*c**3*d**2*e**2*x - 16*b**4*c**3*d*e*(d + e*x)**2
- 8*b**3*c**4*d**4 - 16*b**3*c**4*d**3*e*x + 8*b**3*c**4*d**2*(d + e*x)**2) - 6*c**4*d*e**2*(d + e*x)**(3/2)/(
8*b**7*e**4 - 16*b**6*c*d*e**3 + 16*b**6*c*e**4*x - 48*b**5*c**2*d*e**3*x + 8*b**5*c**2*e**2*(d + e*x)**2 + 16
*b**4*c**3*d**3*e + 48*b**4*c**3*d**2*e**2*x - 16*b**4*c**3*d*e*(d + e*x)**2 - 8*b**3*c**4*d**4 - 16*b**3*c**4
*d**3*e*x + 8*b**3*c**4*d**2*(d + e*x)**2) - 20*c**3*d*e**3*sqrt(d + e*x)/(8*b**6*e**4 - 16*b**5*c*d*e**3 + 16
*b**5*c*e**4*x - 48*b**4*c**2*d*e**3*x + 8*b**4*c**2*e**2*(d + e*x)**2 + 16*b**3*c**3*d**3*e + 48*b**3*c**3*d*
*2*e**2*x - 16*b**3*c**3*d*e*(d + e*x)**2 - 8*b**2*c**4*d**4 - 16*b**2*c**4*d**3*e*x + 8*b**2*c**4*d**2*(d + e
*x)**2) - 6*c**3*d*e*sqrt(d + e*x)/(2*b**6*e**2 - 2*b**5*c*d*e + 2*b**5*c*e**2*x - 2*b**4*c**2*d*e*x) + 6*c**3
*e**3*(d + e*x)**(3/2)/(8*b**6*e**4 - 16*b**5*c*d*e**3 + 16*b**5*c*e**4*x - 48*b**4*c**2*d*e**3*x + 8*b**4*c**
2*e**2*(d + e*x)**2 + 16*b**3*c**3*d**3*e + 48*b**3*c**3*d**2*e**2*x - 16*b**3*c**3*d*e*(d + e*x)**2 - 8*b**2*
c**4*d**4 - 16*b**2*c**4*d**3*e*x + 8*b**2*c**4*d**2*(d + e*x)**2) + 10*c**2*e**4*sqrt(d + e*x)/(8*b**5*e**4 -
16*b**4*c*d*e**3 + 16*b**4*c*e**4*x - 48*b**3*c**2*d*e**3*x + 8*b**3*c**2*e**2*(d + e*x)**2 + 16*b**2*c**3*d*
*3*e + 48*b**2*c**3*d**2*e**2*x - 16*b**2*c**3*d*e*(d + e*x)**2 - 8*b*c**4*d**4 - 16*b*c**4*d**3*e*x + 8*b*c**
4*d**2*(d + e*x)**2) + 4*c**2*e**2*sqrt(d + e*x)/(2*b**5*e**2 - 2*b**4*c*d*e + 2*b**4*c*e**2*x - 2*b**3*c**2*d
*e*x) - 10*d**2*e**2*sqrt(d + e*x)/(-8*b**3*d**4 - 16*b**3*d**3*e*x + 8*b**3*d**2*(d + e*x)**2) + 6*d*e**2*(d
+ e*x)**(3/2)/(-8*b**3*d**4 - 16*b**3*d**3*e*x + 8*b**3*d**2*(d + e*x)**2) - 3*c**2*e**3*sqrt(-1/(c*(b*e - c*d
)**5))*log(-b**3*e**3*sqrt(-1/(c*(b*e - c*d)**5)) + 3*b**2*c*d*e**2*sqrt(-1/(c*(b*e - c*d)**5)) - 3*b*c**2*d**
2*e*sqrt(-1/(c*(b*e - c*d)**5)) + c**3*d**3*sqrt(-1/(c*(b*e - c*d)**5)) + sqrt(d + e*x))/(8*b**2) + 3*c**2*e**
3*sqrt(-1/(c*(b*e - c*d)**5))*log(b**3*e**3*sqrt(-1/(c*(b*e - c*d)**5)) - 3*b**2*c*d*e**2*sqrt(-1/(c*(b*e - c*
d)**5)) + 3*b*c**2*d**2*e*sqrt(-1/(c*(b*e - c*d)**5)) - c**3*d**3*sqrt(-1/(c*(b*e - c*d)**5)) + sqrt(d + e*x))
/(8*b**2) + 3*c**3*d*e**2*sqrt(-1/(c*(b*e - c*d)**5))*log(-b**3*e**3*sqrt(-1/(c*(b*e - c*d)**5)) + 3*b**2*c*d*
e**2*sqrt(-1/(c*(b*e - c*d)**5)) - 3*b*c**2*d**2*e*sqrt(-1/(c*(b*e - c*d)**5)) + c**3*d**3*sqrt(-1/(c*(b*e - c
*d)**5)) + sqrt(d + e*x))/(8*b**3) - 3*c**3*d*e**2*sqrt(-1/(c*(b*e - c*d)**5))*log(b**3*e**3*sqrt(-1/(c*(b*e -
c*d)**5)) - 3*b**2*c*d*e**2*sqrt(-1/(c*(b*e - c*d)**5)) + 3*b*c**2*d**2*e*sqrt(-1/(c*(b*e - c*d)**5)) - c**3*
d**3*sqrt(-1/(c*(b*e - c*d)**5)) + sqrt(d + e*x))/(8*b**3) - c**2*e**2*sqrt(-1/(c*(b*e - c*d)**3))*log(-b**2*e
**2*sqrt(-1/(c*(b*e - c*d)**3)) + 2*b*c*d*e*sqrt(-1/(c*(b*e - c*d)**3)) - c**2*d**2*sqrt(-1/(c*(b*e - c*d)**3)
) + sqrt(d + e*x))/b**3 + c**2*e**2*sqrt(-1/(c*(b*e - c*d)**3))*log(b**2*e**2*sqrt(-1/(c*(b*e - c*d)**3)) - 2*
b*c*d*e*sqrt(-1/(c*(b*e - c*d)**3)) + c**2*d**2*sqrt(-1/(c*(b*e - c*d)**3)) + sqrt(d + e*x))/b**3 + 3*d*e**2*s
qrt(d**(-5))*log(-d**3*sqrt(d**(-5)) + sqrt(d + e*x))/(8*b**3) - 3*d*e**2*sqrt(d**(-5))*log(d**3*sqrt(d**(-5))
+ sqrt(d + e*x))/(8*b**3) - e**2*sqrt(d**(-3))*log(-d**2*sqrt(d**(-3)) + sqrt(d + e*x))/(2*b**3) + e**2*sqrt(
d**(-3))*log(d**2*sqrt(d**(-3)) + sqrt(d + e*x))/(2*b**3) - e*sqrt(d + e*x)/(b**3*d*x) + 3*c**3*d*e*sqrt(-1/(c
*(b*e - c*d)**3))*log(-b**2*e**2*sqrt(-1/(c*(b*e - c*d)**3)) + 2*b*c*d*e*sqrt(-1/(c*(b*e - c*d)**3)) - c**2*d*
*2*sqrt(-1/(c*(b*e - c*d)**3)) + sqrt(d + e*x))/(2*b**4) - 3*c**3*d*e*sqrt(-1/(c*(b*e - c*d)**3))*log(b**2*e**
2*sqrt(-1/(c*(b*e - c*d)**3)) - 2*b*c*d*e*sqrt(-1/(c*(b*e - c*d)**3)) + c**2*d**2*sqrt(-1/(c*(b*e - c*d)**3))
+ sqrt(d + e*x))/(2*b**4) + 3*c*d*e*sqrt(d**(-3))*log(-d**2*sqrt(d**(-3)) + sqrt(d + e*x))/(2*b**4) - 3*c*d*e*
sqrt(d**(-3))*log(d**2*sqrt(d**(-3)) + sqrt(d + e*x))/(2*b**4) + 6*c*e*atan(sqrt(d + e*x)/sqrt(b*e/c - d))/(b*
*4*sqrt(b*e/c - d)) - 6*c*e*atan(sqrt(d + e*x)/sqrt(-d))/(b**4*sqrt(-d)) + 3*c*sqrt(d + e*x)/(b**4*x) - 12*c**
2*d*atan(sqrt(d + e*x)/sqrt(b*e/c - d))/(b**5*sqrt(b*e/c - d)) + 12*c**2*d*atan(sqrt(d + e*x)/sqrt(-d))/(b**5*
sqrt(-d))

________________________________________________________________________________________

Giac [B] Leaf count of result is larger than twice the leaf count of optimal. 508 vs. \(2 (226) = 452\).
time = 1.45, size = 508, normalized size = 2.07 \begin {gather*} -\frac {{\left (48 \, c^{4} d^{2} - 84 \, b c^{3} d e + 35 \, b^{2} c^{2} e^{2}\right )} \arctan \left (\frac {\sqrt {x e + d} c}{\sqrt {-c^{2} d + b c e}}\right )}{4 \, {\left (b^{5} c d - b^{6} e\right )} \sqrt {-c^{2} d + b c e}} + \frac {24 \, {\left (x e + d\right )}^{\frac {7}{2}} c^{4} d^{2} e - 72 \, {\left (x e + d\right )}^{\frac {5}{2}} c^{4} d^{3} e + 72 \, {\left (x e + d\right )}^{\frac {3}{2}} c^{4} d^{4} e - 24 \, \sqrt {x e + d} c^{4} d^{5} e - 24 \, {\left (x e + d\right )}^{\frac {7}{2}} b c^{3} d e^{2} + 108 \, {\left (x e + d\right )}^{\frac {5}{2}} b c^{3} d^{2} e^{2} - 144 \, {\left (x e + d\right )}^{\frac {3}{2}} b c^{3} d^{3} e^{2} + 60 \, \sqrt {x e + d} b c^{3} d^{4} e^{2} + {\left (x e + d\right )}^{\frac {7}{2}} b^{2} c^{2} e^{3} - 40 \, {\left (x e + d\right )}^{\frac {5}{2}} b^{2} c^{2} d e^{3} + 85 \, {\left (x e + d\right )}^{\frac {3}{2}} b^{2} c^{2} d^{2} e^{3} - 46 \, \sqrt {x e + d} b^{2} c^{2} d^{3} e^{3} + 2 \, {\left (x e + d\right )}^{\frac {5}{2}} b^{3} c e^{4} - 13 \, {\left (x e + d\right )}^{\frac {3}{2}} b^{3} c d e^{4} + 9 \, \sqrt {x e + d} b^{3} c d^{2} e^{4} + {\left (x e + d\right )}^{\frac {3}{2}} b^{4} e^{5} + \sqrt {x e + d} b^{4} d e^{5}}{4 \, {\left (b^{4} c d^{2} - b^{5} d e\right )} {\left ({\left (x e + d\right )}^{2} c - 2 \, {\left (x e + d\right )} c d + c d^{2} + {\left (x e + d\right )} b e - b d e\right )}^{2}} + \frac {{\left (48 \, c^{2} d^{2} - 12 \, b c d e - b^{2} e^{2}\right )} \arctan \left (\frac {\sqrt {x e + d}}{\sqrt {-d}}\right )}{4 \, b^{5} \sqrt {-d} d} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((e*x+d)^(1/2)/(c*x^2+b*x)^3,x, algorithm="giac")

[Out]

-1/4*(48*c^4*d^2 - 84*b*c^3*d*e + 35*b^2*c^2*e^2)*arctan(sqrt(x*e + d)*c/sqrt(-c^2*d + b*c*e))/((b^5*c*d - b^6
*e)*sqrt(-c^2*d + b*c*e)) + 1/4*(24*(x*e + d)^(7/2)*c^4*d^2*e - 72*(x*e + d)^(5/2)*c^4*d^3*e + 72*(x*e + d)^(3
/2)*c^4*d^4*e - 24*sqrt(x*e + d)*c^4*d^5*e - 24*(x*e + d)^(7/2)*b*c^3*d*e^2 + 108*(x*e + d)^(5/2)*b*c^3*d^2*e^
2 - 144*(x*e + d)^(3/2)*b*c^3*d^3*e^2 + 60*sqrt(x*e + d)*b*c^3*d^4*e^2 + (x*e + d)^(7/2)*b^2*c^2*e^3 - 40*(x*e
+ d)^(5/2)*b^2*c^2*d*e^3 + 85*(x*e + d)^(3/2)*b^2*c^2*d^2*e^3 - 46*sqrt(x*e + d)*b^2*c^2*d^3*e^3 + 2*(x*e + d
)^(5/2)*b^3*c*e^4 - 13*(x*e + d)^(3/2)*b^3*c*d*e^4 + 9*sqrt(x*e + d)*b^3*c*d^2*e^4 + (x*e + d)^(3/2)*b^4*e^5 +
sqrt(x*e + d)*b^4*d*e^5)/((b^4*c*d^2 - b^5*d*e)*((x*e + d)^2*c - 2*(x*e + d)*c*d + c*d^2 + (x*e + d)*b*e - b*
d*e)^2) + 1/4*(48*c^2*d^2 - 12*b*c*d*e - b^2*e^2)*arctan(sqrt(x*e + d)/sqrt(-d))/(b^5*sqrt(-d)*d)

________________________________________________________________________________________

Mupad [B]
time = 2.01, size = 2500, normalized size = 10.20 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((d + e*x)^(1/2)/(b*x + c*x^2)^3,x)

[Out]

(((d + e*x)^(3/2)*(b^4*e^5 + 72*c^4*d^4*e - 144*b*c^3*d^3*e^2 + 85*b^2*c^2*d^2*e^3 - 13*b^3*c*d*e^4))/(4*b^4*(
c*d^2 - b*d*e)) - ((d + e*x)^(1/2)*(b^3*e^4 + 24*c^3*d^3*e - 36*b*c^2*d^2*e^2 + 10*b^2*c*d*e^3))/(4*b^4) + ((b
*e - 2*c*d)*(d + e*x)^(5/2)*(b^2*c*e^3 + 18*c^3*d^2*e - 18*b*c^2*d*e^2))/(2*b^4*(c*d^2 - b*d*e)) + (c*e*(d + e
*x)^(7/2)*(24*c^3*d^2 + b^2*c*e^2 - 24*b*c^2*d*e))/(4*b^4*(c*d^2 - b*d*e)))/(c^2*(d + e*x)^4 - (d + e*x)*(4*c^
2*d^3 + 2*b^2*d*e^2 - 6*b*c*d^2*e) - (4*c^2*d - 2*b*c*e)*(d + e*x)^3 + (d + e*x)^2*(b^2*e^2 + 6*c^2*d^2 - 6*b*
c*d*e) + c^2*d^4 + b^2*d^2*e^2 - 2*b*c*d^3*e) - (atan((((((d + e*x)^(1/2)*(b^6*c^3*e^8 + 4608*c^9*d^6*e^2 - 13
824*b*c^8*d^5*e^3 + 22*b^5*c^4*d*e^7 + 15072*b^2*c^7*d^4*e^4 - 7104*b^3*c^6*d^3*e^5 + 1226*b^4*c^5*d^2*e^6))/(
8*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)) + (((b^14*c^2*d*e^7 - 24*b^10*c^6*d^5*e^3 + 60*b^11*c^5*d^4*e^
4 - 46*b^12*c^4*d^3*e^5 + 9*b^13*c^3*d^2*e^6)/(b^12*c^2*d^4 + b^14*d^2*e^2 - 2*b^13*c*d^3*e) - ((d + e*x)^(1/2
)*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e)*(128*b^10*c^5*d^5*e^2 - 320*b^11*c^4*d^4*e^3 + 256*b^12*c^3*d^3*e^4 - 64
*b^13*c^2*d^2*e^5))/(64*b^5*(d^3)^(1/2)*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)))*(b^2*e^2 - 48*c^2*d^2 +
12*b*c*d*e))/(8*b^5*(d^3)^(1/2)))*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e)*1i)/(8*b^5*(d^3)^(1/2)) + ((((d + e*x)^
(1/2)*(b^6*c^3*e^8 + 4608*c^9*d^6*e^2 - 13824*b*c^8*d^5*e^3 + 22*b^5*c^4*d*e^7 + 15072*b^2*c^7*d^4*e^4 - 7104*
b^3*c^6*d^3*e^5 + 1226*b^4*c^5*d^2*e^6))/(8*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)) - (((b^14*c^2*d*e^7
- 24*b^10*c^6*d^5*e^3 + 60*b^11*c^5*d^4*e^4 - 46*b^12*c^4*d^3*e^5 + 9*b^13*c^3*d^2*e^6)/(b^12*c^2*d^4 + b^14*d
^2*e^2 - 2*b^13*c*d^3*e) + ((d + e*x)^(1/2)*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e)*(128*b^10*c^5*d^5*e^2 - 320*b^
11*c^4*d^4*e^3 + 256*b^12*c^3*d^3*e^4 - 64*b^13*c^2*d^2*e^5))/(64*b^5*(d^3)^(1/2)*(b^8*c^2*d^4 + b^10*d^2*e^2
- 2*b^9*c*d^3*e)))*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e))/(8*b^5*(d^3)^(1/2)))*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*
e)*1i)/(8*b^5*(d^3)^(1/2)))/((1728*c^10*d^6*e^3 - (35*b^6*c^4*e^9)/32 - 5184*b*c^9*d^5*e^4 + (63*b^5*c^5*d*e^8
)/4 + 5508*b^2*c^8*d^4*e^5 - 2376*b^3*c^7*d^3*e^6 + (1233*b^4*c^6*d^2*e^7)/4)/(b^12*c^2*d^4 + b^14*d^2*e^2 - 2
*b^13*c*d^3*e) + ((((d + e*x)^(1/2)*(b^6*c^3*e^8 + 4608*c^9*d^6*e^2 - 13824*b*c^8*d^5*e^3 + 22*b^5*c^4*d*e^7 +
15072*b^2*c^7*d^4*e^4 - 7104*b^3*c^6*d^3*e^5 + 1226*b^4*c^5*d^2*e^6))/(8*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*
c*d^3*e)) + (((b^14*c^2*d*e^7 - 24*b^10*c^6*d^5*e^3 + 60*b^11*c^5*d^4*e^4 - 46*b^12*c^4*d^3*e^5 + 9*b^13*c^3*d
^2*e^6)/(b^12*c^2*d^4 + b^14*d^2*e^2 - 2*b^13*c*d^3*e) - ((d + e*x)^(1/2)*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e)*
(128*b^10*c^5*d^5*e^2 - 320*b^11*c^4*d^4*e^3 + 256*b^12*c^3*d^3*e^4 - 64*b^13*c^2*d^2*e^5))/(64*b^5*(d^3)^(1/2
)*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)))*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e))/(8*b^5*(d^3)^(1/2)))*(b^
2*e^2 - 48*c^2*d^2 + 12*b*c*d*e))/(8*b^5*(d^3)^(1/2)) - ((((d + e*x)^(1/2)*(b^6*c^3*e^8 + 4608*c^9*d^6*e^2 - 1
3824*b*c^8*d^5*e^3 + 22*b^5*c^4*d*e^7 + 15072*b^2*c^7*d^4*e^4 - 7104*b^3*c^6*d^3*e^5 + 1226*b^4*c^5*d^2*e^6))/
(8*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)) - (((b^14*c^2*d*e^7 - 24*b^10*c^6*d^5*e^3 + 60*b^11*c^5*d^4*e
^4 - 46*b^12*c^4*d^3*e^5 + 9*b^13*c^3*d^2*e^6)/(b^12*c^2*d^4 + b^14*d^2*e^2 - 2*b^13*c*d^3*e) + ((d + e*x)^(1/
2)*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e)*(128*b^10*c^5*d^5*e^2 - 320*b^11*c^4*d^4*e^3 + 256*b^12*c^3*d^3*e^4 - 6
4*b^13*c^2*d^2*e^5))/(64*b^5*(d^3)^(1/2)*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)))*(b^2*e^2 - 48*c^2*d^2
+ 12*b*c*d*e))/(8*b^5*(d^3)^(1/2)))*(b^2*e^2 - 48*c^2*d^2 + 12*b*c*d*e))/(8*b^5*(d^3)^(1/2))))*(b^2*e^2 - 48*c
^2*d^2 + 12*b*c*d*e)*1i)/(4*b^5*(d^3)^(1/2)) - (atan((((-c^3*(b*e - c*d)^3)^(1/2)*(((d + e*x)^(1/2)*(b^6*c^3*e
^8 + 4608*c^9*d^6*e^2 - 13824*b*c^8*d^5*e^3 + 22*b^5*c^4*d*e^7 + 15072*b^2*c^7*d^4*e^4 - 7104*b^3*c^6*d^3*e^5
+ 1226*b^4*c^5*d^2*e^6))/(8*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)) + ((-c^3*(b*e - c*d)^3)^(1/2)*((b^14
*c^2*d*e^7 - 24*b^10*c^6*d^5*e^3 + 60*b^11*c^5*d^4*e^4 - 46*b^12*c^4*d^3*e^5 + 9*b^13*c^3*d^2*e^6)/(b^12*c^2*d
^4 + b^14*d^2*e^2 - 2*b^13*c*d^3*e) - ((-c^3*(b*e - c*d)^3)^(1/2)*(d + e*x)^(1/2)*(35*b^2*e^2 + 48*c^2*d^2 - 8
4*b*c*d*e)*(128*b^10*c^5*d^5*e^2 - 320*b^11*c^4*d^4*e^3 + 256*b^12*c^3*d^3*e^4 - 64*b^13*c^2*d^2*e^5))/(64*(b^
8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)*(b^8*e^3 - b^5*c^3*d^3 + 3*b^6*c^2*d^2*e - 3*b^7*c*d*e^2)))*(35*b^2*
e^2 + 48*c^2*d^2 - 84*b*c*d*e))/(8*(b^8*e^3 - b^5*c^3*d^3 + 3*b^6*c^2*d^2*e - 3*b^7*c*d*e^2)))*(35*b^2*e^2 + 4
8*c^2*d^2 - 84*b*c*d*e)*1i)/(8*(b^8*e^3 - b^5*c^3*d^3 + 3*b^6*c^2*d^2*e - 3*b^7*c*d*e^2)) + ((-c^3*(b*e - c*d)
^3)^(1/2)*(((d + e*x)^(1/2)*(b^6*c^3*e^8 + 4608*c^9*d^6*e^2 - 13824*b*c^8*d^5*e^3 + 22*b^5*c^4*d*e^7 + 15072*b
^2*c^7*d^4*e^4 - 7104*b^3*c^6*d^3*e^5 + 1226*b^4*c^5*d^2*e^6))/(8*(b^8*c^2*d^4 + b^10*d^2*e^2 - 2*b^9*c*d^3*e)
) - ((-c^3*(b*e - c*d)^3)^(1/2)*((b^14*c^2*d*e^7 - 24*b^10*c^6*d^5*e^3 + 60*b^11*c^5*d^4*e^4 - 46*b^12*c^4*d^3
*e^5 + 9*b^13*c^3*d^2*e^6)/(b^12*c^2*d^4 + b^14...

________________________________________________________________________________________